Primordial nuggets survival and QCD pairing 
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We revisit the problem of boiling and surface evaporation of quark nuggets in the cosmological 
quark-hadron transition with the explicit consideration of pairing between quarks in a color-flavor 
locked (CFL) state. Assuming that primordial quark nuggets are actually formed, we analyze the 
consequences of pairing on the rates of boiling and surface evaporation in order to determine whether 
they could have survived with substantial mass. We find a substantial quenching of the evaporation 
+ boiling processes, which suggests the survival of primordial nuggets for the currently considered 
range of the pairing gap A. Boiling is shown to depend on the competition of an increased stability 
window and the suppression of the rate, and is not likely to dominate the destruction of the nuggets. 
If surface evaporation dominates, the fate of the nuggets depend on the features of the initial mass 
spectrum of the nuggets, their evaporation rate, and the value of the pairing gap, as shown and 
discussed in the text. 



I. INTRODUCTION 

Approximately 10 -5 seconds after the Big Bang the 
early universe was filled with a hot and expanding mix- 
ture of elementary particles. The Universe was composed 
mainly by photons, charged leptons, neutrinos, quarks 
and gluons (and the corresponding antiparticles) coex- 
isting in thermal and chemical equilibrium through elec- 
troweak interactions. As the Universe expanded, this 
mixture cooled down to a critical temperature at which 
the plasma of free quarks and gluons converted into 
hadrons. Early studies of this transition started in the 
1980's [Hi and gave a broad-brush picture of the 
physics involved (for a more complete reference list see 

EliHii). 

A very important question is whether the transition 
is actually first order, second order or just a crossover. 
Dramatic effects are expected in the first case, while a 
second order or crossover would be much less spectacu- 
lar. Lattice numerical simulation is the best approach 
currently available for the study of QCD near the fi- 
nite temperature transition point. While it has longly 
been known that the transition is first order in the case 
of pure gluonic calculations (corresponding to infinitely 
heavy quarks), and in the case of four light quarks, the 
actual physical case is elusive. At present, there are well 
established non-perturbative lattice techniques to study 
this transition at fi = and T^O. The order of the tran- 
sition is known as a function of the quark masses showing 
that the physical point is probably in the crossover re- 
gion unless the s quark mass is small (in which case it 
should be first order). For recent reviews see 0,^3 an d 
references therein. 

Interesting baryon fluctuations would have been pro- 
duced by a first order transition. The two phases need to 
coexist long enough for baryon transport to shuffle the 
baryon number across the phase boundary. As pointed 
out in early studies, the onset of the supposedly first- 
order transition requires some degree of supercooling Q . 



If the transition is not first order, no supercooling could 
possibly occur (even if the equation of state gave rise to 
a very rapid change in the energy density) due to the 
extremely slow expansion of the Universe. 

The generation of primordial isothermal baryon num- 
ber inhomogcncitics can be understood within a scenario 
of cosmic separation of phases 0, || • When the uni- 
verse cools to the critical temperature Tqcd nucleation 
of bubbles of the hadron phase could begin. However, 
it is a general feature of the nucleation theory that the 
nucleation probability is not large enough at the criti- 
cal temperature but for temperatures below it. There- 
fore, the universe supercools below Tqcd being still in 
the quark phase until the nucleation rate becomes suf- 
ficiently large. After a brief stage of nucleation during 
which the hadron bubbles grow and reheat the universe 
back to Tqcd, nucleation is again inhibited due to its 
low probability and the expansion of the universe makes 
hadron bubbles to grow slowly at expenses of the quark 
phase. Once hadron bubbles occupy roughly half of the 
total volume they are able to collide and merge leav- 
ing the universe with shrinking droplets of quark-gluon 
plasma immersed in a hadron matter medium. 

The fate of these baryon number inhomogeneities de- 
pend on how heat and baryon number are transported 
across the transition front |j, |8|. Latent heat (or en- 
tropy) could be carried out by neutrinos, surface evap- 
oration of hadrons (mostly pions) and by the motion of 
the transition front which converts volume of one vacuum 
into another. The baryon number transport across the 
conversion front depends on the bulk properties of both 
phases and on the penetrability of the interface (which 
quantifies the chance of baryon number to pass from one 
side of the boundary to the other) . Estimations of baryon 
number penetrability have been made within the frame of 
the chromoelcctric flux tube model 
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baryon penetrability indeed happens to be small, it may 
be possible to accumulate almost all the baryon number 
density in the quark-gluon phase (see below). In such 
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a case, and depending on the parameters, the inhomo- 
geneities may be large enough to produce strange quark 
matter (SQM). This results in a universe in thermal equi- 
librium but with an inhomogcncous baryon distribution 
(i.e. out of chemical equilibrium). 

The study of the extreme case in which quark nug gets 
form has been undertaken by a number of authors [L3 . 
EmEHIHIIllEl. An absolute upper limit to the 
baryon number contained in the nuggets is determined 
by the size of the cosmological horizon evaluated at the 
critical temperature. The simplest estimate yields the 
well-known value 
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Actually, the details of the dynamics will determine an 
initial mass function at the end of the transition. This 
is a quite complicated problem and has not been solved 
in detail, although Bhattacharyya et al. [23| presented a 
series of calculations showing that the maximum baryon 
number of the nuggets is ~ 10 43 for Tqcd — 150 MeV, 
which fit comfortably within the horizon size. 

After the QCD phase transition, the temperature 
in the primordial Universe is still high enough to al- 
low for evaporation of hadrons from the surface of the 
nuggets, and in principle to allow for boiling (nucleation 
of hadronic bubbles inside its volume). This is a con- 
sequence of the presence of the —TS term in the free 
energy, which disfavors the strange quark matter phase 
at intermediate temperatures. It is only at low temper- 
atures (T w 2 MeV) that nuggets begin to be preferred 
to free hadrons. Previous work found that boiling is not 
possible for reasonable values of the bag constant B since 
the timescale is too short for bubble nucleation to take 
place (23, [2l| . If so, surface evaporation seems to be the 
only mechanism able to destroy the primordial nuggets, 
although the very survival of these entities may be con- 
sidered as still subject to uncertainties. Since it is likely 
that quarks inside the nuggets settle in paired states at a 
relatively high temperature (see Fig. 1 for a qualitative 
sketch), we shall examine in the remaining of this work 
the effects of quark pairing on the evaporation/boiling at 
intermediates temperatures, thus revisiting the question 
of nugget survival. 



II. BOILING OF CFL NUGGETS 

As stated above, quark nuggets are born hot and there- 
fore nucleation of hadronic bubbles could occur inside 
them. Nevertheless, as it has been shown previously [2lj| . 
boiling is unlikely to destroy primordial nuggets for rea- 
sonable values of the bag constant B. However, we must 
note that pairing must occur at temperatures below the 
critical temperature Xa = 0.57A. Therefore, the anal- 
ysis made in |2l| holds only in the temperature regime 
between Tqcd and Ta while below Ta pairing effects 



> 

CD 



-150 



QGP 








A 




\b 






Hadrons \^ 


( 


CFL 


. preferred 





-300 



|i [MeV] 



FIG. 1: Path of the strange quark matter nuggets in the T — fj, 
plane. Nuggets are quickly formed starting at T = Tqcd ~ 
150 MeV, and they are fragile to evaporation/boiling at in- 
termediate temperatures as discussed in text. A transition 
to the CFL phase occurs at the points marked with crosses. 
The path labelled as "A" assumes a very quick formation 
of the nugget (that is, t formation < H' 1 , see Ref. [23]), 
which evolves at constant density afterwards following a ver- 
tical path. A (perhaps more realistic) path "B" has been also 
sketched, in which the formation is slower tf orma uon ~ H -1 . 
After crossing the CFL boundary nuggets are "safe" because 
the pairing now protects them against evaporation/boiling 
and attaining the dashed line temperature is no longer rele- 
vant for their fate. Thus, their masses freeze at a higher value 
when quarks become locked in CFL states. 



must be taken into account. A remarkable consequence 
of QCD pairing is that the stability window for strange 
matter is considerably enlarged, allowing a wider range 
for B 0. Thus, although pairing should difficult boiling 
because more energy is necessary to produce an hadron 
lump, it is not clear a priori to what extent the modifi- 
cation of the stability characteristics of SQM counteracts 
this effect. 

Let us briefly examine the thermodynamic description 
of boiling of primordial nuggets including the effect of 
pairing between quarks. We assume, for simplicity, that 
the nucleated phase is in thermal and chemical equilib- 
rium with the nugget, which itself evolves at fixed /1 (that 
is, along a path of the type "A" in Fig. 1). The work done 
to form a bubble of radius r composed by the hadronic 
phase inside the quark phase is 



W = --nr 3 AP + 47TOT 2 - 2njr + -irr 3 x -n B A, (2) 
3 3 2 

being AP = Ph — P q the pressure difference between both 
phases, a = a q + <Jh the surface tension, 7 = j q — 7^ the 
curvature coefficient, jig the baryon number density in 
the hadronic phase, and A the gap of the CFL pair. The 
innovation here is the last term which is introduced by 
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considering that to put three quarks in an hadron an en- 
ergy A must be expended for each CFL pair that is dis- 
assembled inside the quark phase. Further refinements 
should be considered if hadrons are assumed to be com- 
posed by a diquark plus a free quark; in this case the 
net energy released (per baryon) should be ~ ^A. The 
effect of A on a and 7 themselves is unknown and will be 
neglected in this first approach. Note that the gap also 
enters the free energy trough the pressure as a quadratic 
term. The critical radius r c is found by extremizing W, 



being b = "/F/2a 2 and 



F = AP - -n B A. 
2 



(3) 



(4) 



Only those bubbles with a radius greater than r c will 
be able to grow. The nucleation rate for the growing 
bubbles is given by 



K boU =T 4 exp(-W c /T) 



(5) 



being W c = ^-[2 + 2(1 - bf 2 - 36]. As we shall see 
below the contribution of A tends to suppress the rate 
since it enters in such a way that F becomes smaller and 
W c becomes larger. However, since the stability behavior 
is modified by pairing, allowing stability for a much wider 
rage of the bag constant P, it is necessary to determine 
which is the leading effect on the boiling process, as we 
shall do in the following. 

The effect of boiling on the quark nuggets can be ana- 
lyzed by means of a slightly different condition (see pl| 
and references therein). Comparing the total area of the 
nugget and the bubbles 



From an inspection of Eq. (JJJ it is easily recognized 
that the main effect of pairing on boiling enters only 
through a boost in the bag constant B. This can be un- 
derstood by comparing the boiling of CFL strange matter 
with the boiling of unpaired SQM. If we assume that the 
strange quark mass m s is zero, the pressure of SQM and 
CFL strange matter differ only by a term proportional 
to A 2 , i.e. P S qm = Pcfl - (3AV)A 2 M- In both, 
the CFL and the unpaired SQM phases, there exists a 
symmetric flavor composition with n u = nd — n s . Note 
that the last will be not true when considering that the 
strange quark has a finite mass; in this case the CFL 
phase will still have a symmetric composition but SQM 
will not. However, we do not expect strong departures 
from this simple analysis. The nucleation of hadron bub- 
bles occurs in a very fast timescale r s ~ 10~ 24 s, typical of 
strong interactions, therefore the transition happens out 
of weak equilibrium irrespective of the pairing of quark 
matter. This means that flavor will be conserved during 
the process of nucleation (only after t w ~ 1CP 8 s, will 
/3-decays lead the just formed hadron phase to its equi- 
librium configuration) . The conservation of flavor during 
the formation of bubbles, and the fact of both SQM and 
CFL strange matter having the same composition (in a 
first approximation) guarantees that exactly the same gas 
of hadrons will be produced both beginning with a CFL 
or with an unpaired SQM composition. Therefore, the 
pressure difference Ap) between the SQM phase and the 
hadron phase, and the difference AP between the CFL 
phase and the hadron one, are related by 



AP = APn 



3AV 



2, ,2 



(8) 



The simple relation given by Eq. JSJl allow us to gauge 
straightforwardly the impact of pairing in the boiling pro- 
cess. The difference F can be written, in the case of 
massless quarks, as 
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it is found that boiling is less important than surface 
evaporation for a baryon number A below the value Ab u 
given by 



F = AP - -n B A 
2 

= Ph. — Pfree — B, 



APo - -iAV ~ ~n B A 



} eff- 



(9) 



The effect of pairing has been included in the bag con- 
stant B by defining an "effective bag constant" 
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; [2 + 2(1 - fo) 3 / 2 - 36] ; 



TF 2 



( 7 ) 



For a given value of AboU, the last equation gives the 
critical bag constant B and er separating the boiling and 
the non-boiling regions. The value of the critical B and 
a is almost insensitive to the value of A^ou . Spanning the 
range 1 < Ab u < 10 57 only c hang es the critical values of 
P 1 / 4 and a 1 / 3 by a few MeV [lj 



B, 



B - ^AV 



(10) 



where Pfree is the pressure of a flavor-symmetric mixture 
of free quarks. Note that pairing enters trough a lead- 
ing contribution ~ fj, 3 A associated with the condensation 
work and a second order contribution ~ A 2 /i 2 associated 
with the modification of the pressure in the CFL phase. 
Although our analysis does not include the important ef- 
fects of the finite mass of the strange quark and of the 
chemical composition of the phases, it seems clear that 
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the here shown tendency should be qualitatively the same 
in a more realistic study. Also, as stated above, some re- 
finements would need to be considered due to the effect 
of finite temperature, and the pairing gap in the sur- 
face tension and curvature terms. In summary, all the 
difference with the boiling of quark nuggets made up of 
unpaired massless quarks is that we must use here the 
effective bag constant defined by Eq. (JTDJ. 

The likelihood of boiling can be analyzed in the B 
vs. T plane (see Fig 2). For an unpaired quark mix- 
ture the critical B above which boiling is allowed (for 
any baryon number A) is always greater than the max- 
imum B that allows stable SQM for a transition out of 
the equilibrium [2]]]. As we have shown in Eq. (|10|l . 
the pairing shifts up the critical curve by a magnitude 
(3A 2 /x 2 )/7r 2 + (3n B A)/2. On the other hand, as shown 
in Ref. |24|, the maximum B that allows stable CFL 
strange matter also shifts up a magnitude m 2 A 2 /(37r 2 ). 
Therefore, the net shift h& is 



h\ = 



(V~m 2 )A 2 
3tt 2 



3 A 

2 nsA ' 



(11) 



which is clearly positive in the range of interest since the 
leading term is the second one and scales as /i 3 A. This 
means that nucleation is impossible in the temperature 
regime where QCD pairing operates (T < Ta), even in 
the most favorable situation in which the nucleated phase 
is in equilibrium. 
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FIG. 2: Sketch of the effects of CFL pairing on the boiling 
of nuggets. The stability window of SQM is realized between 
the lower dotted line and the solid horizontal line. Analo- 
gously, the wider stability window of CFL SQM holds between 
the lower dotted line and the short-dashed horizontal line as 
marked |24j. For a given temperature, boiling of primordial 
quark nuggets is allowed only above the parabolic-like curves 
(solid line for unpaired SQM and dashed for CFL SQM). In 
spite of the rising of the upper stability limit for CFL SQM, 
the boiling curve also raises, and since both curves lie above 
the corresponding upper limits for matter stability (horizon- 
tal lines) , quark matter nuggets must survive boiling (with or 
without pairing). 



III. SURFACE EVAPORATION 

Surface evaporation of hadrons has been addressed as 
a first mechanism for nugget destruction by Alcock and 
Farhi 0] as explained above. Using simple detailed bal- 
ance arguments, they concluded that the flux of baryons 
from the surface was more than enough to evaporate the 
nuggets for all but the highest (unphysical) masses. Fur- 
ther work revisited the issue by employing chromoelectric 
flux tube expressions |l3j . which happened to be much 
smaller than the naive flux employed originally. It was 
found that nuggets having a baryon number larger than 
10 39 could survive evaporation. Another detailed study 
by Madsen, Heiselberg and Riisager also considered 
explicitly the effect of flavor non-equilibrium at the sur- 
face of the nugget and found that lumps with baryon 
number as low as A ~ 10 46 could survive evaporation. 
These and other calculations suggests that a low sur- 
face baryon flux allows the survival of large, but not ex- 
treme nuggets, perhaps down to A ~ 10 40 provided the 
evaporation flux is low enough. Therefore it is of inter- 
est to understand what happens when quark pairing is 
introduced in this picture. 

An evaporating lump is slightly cooler than the en- 
vironment, which for temperatures ~ 100 MeV is com- 
posed mainly by photons, neutrinos, electrons and their 
antiparticles. Heat flows from the surrounding medium 



into the lump, providing the energy to power the evap- 
oration. Near the surface of the lump there is a high 
concentration of baryons that have just evaporated. Due 
to their large mean free path, neutrinos are the most 
efficient energy carriers. We shall discuss the effects of 
pairing using a scaling of the simplest rate derived by 
Alcock and Farhi |l5|. since the latter provides a good 
description of all surface evaporation rates presented in 
the literature. The number of hadrons evaporated from 
the surface per second is written as 

TZ = a ^Uhl T 2 A 2/3 e ~i/T (12) 

with I = 20 MeV the binding energy, a = 3.1 x 10~ 4 
MeV -2 and m„ the neutron mass. The parameter a 
is introduced in order to reproduce approximately the 
behavior of the flux within very different models, which 
differ by several orders of magnitude (see e.g. 0, Q, 

E3). 

While above the critical temperature for pair forma- 
tion Ta, the evaporation rate would be given by Ea. (|12ll . 
Below Ta, we use the same combinatorial criterion as in 
the previous section for the breakup of quark pairs to 
write down the rate of evaporation from CFL nuggets 
7?.A in a first approximation as 
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FIG. 3: The baryon number evolution of quark nuggets as a 
function of the temperature of the universe for A = 50 MeV. 
The upper dashed line is the baryon number included within 
the horizon. The lower dashed line (given by G 2 F T 4 A 1/3 » 1) 
divides the regions of diffusive neutrino heating and neutrino 
transparency of the environment of the nuggets. There is al- 
most no evaporation in the diffusive neutrino heating regime, 
simply because not enough energy is supplied to power the 
evaporation. The vertical line shows the temperature Ta be- 
low which pairing operates. The solid lines show the evolu- 
tion of nuggets with three different initial baryon numbers. 
The heavier nugget reaches Ta even before leaving the diffu- 
sive neutrino heating regime, and therefore freezes out retain- 
ing its initial baryon number. The lighter of these nuggets 
evaporates completely as soon as it enters the neutrino trans- 
parent regime. The intermediate mass nugget maintains its 
initial baryon number until it enters the neutrino transpar- 
ent regime. Thereafter, it evaporates substantially until it 
reaches the critical temperature Ta where it freezes out with 
a smaller mass. This nugget would have evaporated in the 
absence of pairing, as is apparent from the dotted curve cor- 
responding to A = 0. These calculations assume a — 10~ 3 , 
which is much higher than the values given by chromoelectric 
flux tube models but smaller than the extreme value given by 
detailed balance. 



TZ A = K X e"3A/2T 



(13) 



Since the energy cost of pair breakup is a general fea- 
ture of the models, we expect this suppression to hold 
quite independently of the detailed physics. The impor- 
tant feature is that surface evaporation rates get effec- 
tively quenched as soon as the CFL phase sets in, at a 
temperature Ta = 0.57A, which is certainly much higher 
than the ~ 2 MeV necessary to stabilize the nuggets. 
Therefore, it may be said that CFL states freeze out the 
mass of the nuggets once they cool down to Ta . 

The parametrization of the surface evaporation rates 
given in Eqs. Ijl2(l and 1|13|1 allows a simple analytic solu- 
tion of the evolution equation dA/dt = TZ of the baryon 
number of the nugget. This solution is 
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FIG. 4: The same as the previous figure but for A = 100 
MeV. Lower masses can now reach the pairing temperature 
and freeze out 



A l / 3 (T) = A, 
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where we have identified the initial baryon number Aq = 
A(T ) and C = a(-2a m n ) / (6tt 2 ) x [45/(172tt 3 G')] 1 / 2 
(with G the Newton constant). 

As discussed in Refs. 0, the beginning of the 
evaporation is possible when the nuggets are surrounded 
by an optically thin environment. The neutrino influx 
is then capable of powering the baryon evaporation at 
the nugget surface. The transition from the diffusive to 
transparent regime satisfies the condition GpT^A 1 / 3 w 
1 (with Gf the Fermi constant). Evaporation is small 
before the nugget crosses this boundary curve. Thus, the 
temperature To at which each nugget begin to evaporate 
is a function of the initial baryon content Ao, assumed 
to be the value it had at the formation temperature (see 
Figs. 3 and 4). 

It can be checked that, depending on the value of the 
pairing gap A, a subset of nuggets that proceed directly 
from the diffusive neutrino heating regime to the CFL 
paired state may exist. Therefore, these nuggets remain 
essentially frozen with the same baryon number they 
had at their formation. The minimum baryon number 
A freeze that satisfies this condition is found inserting the 
relation Ta = 0.57A in the transition condition above, 
and is given by 
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Nuggets smaller than Af reeze will evaporate substan- 
tially once they enter the optically thin neutrino regime, 
but may survive if they manage to cool down to Ta with 
some finite mass. 
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While quarks remain unpaired, the evaporation rate 
will be given by Eq. fHft. Therefore, for T A < T < T 
the baryon number density as a function of temperature 
is given by 



A 1 / 3 (T) = Al /3 - C [Ei(-I/T) - Ei{-I/T )] (16) 

where Ei{x) is the exponential integral. 

Once T < Ta, pairing reduces the evaporation rate to 
the expression Eg . (|13|) . and the baryon number density 
as a function of temperature follows 



A 1/3 {T) = Al /3 



C 



Ei 



-I 



- Ei 



-I 



T J ~ \ T A 
Ei(-I/T A )-Ei(-I/T Q ) . 



(17) 



A simple approximation for Ei (x) which is good within 
a few percent in the range of interest, is the following 



Ei{x) 



exp(ir) 
x 2 — 2x 



(18) 



which is useful for understanding the relative weight of 
each term in the corresponding temperature regimes. 
The complete results are depicted in Figs. 3 and 4. The 
effects of CFL pairing are apparent when nuggets reach 
Ta, since many of them are able to survive while they 
would have been evaporated in the absence of this pair- 
ing. As a corollary, we may state quite generally that a 
given initial mass function of nuggets would be stretched 
towards the smallest masses because of CFL pairing. De- 
tailed calculations of this features will be presented in a 
future publication. We finally stress that this evaporat- 
ing population may not exist at all, depending on the 
form of the initial mass function. 



IV. CONCLUSIONS 



may survive if their destruction is not complete when the 
Universe cools down to a sufficiently low temperature. 

We have shown in this work that the consideration 
of pairing brings an additional twist to the problem of 
nugget survival at intermediate temperatures. Specifi- 
cally, we have shown that both the boiling and the sur- 
face evaporation get suppressed because of the presence 
of the gap A in the respective rates. 

Boiling of nuggets has been already discussed in the 
literature and found unlikely in the most realistic calcu- 
lations. When CFL pairing is included, the boiling is also 
unlikely because, in spite of the increase of the stability 
window, the rate is suppressed by A and the net effect 
produces h& > in realistic cases. 

In the case of surface evaporation, the fate of the 
nuggets depend mainly on the (unknown) characteristics 
of the initial mass spectrum of the nuggets, their evapo- 
ration rate, and the value of the pairing gap. However, 
and independently of these uncertainties, many general 
trends can be noticed. If the value of the pairing gap 
A is sufficiently high, the nuggets perhaps as small as 
~ 10 42 and up to A max enter the CFL phase before leav- 
ing the regime which is opaque to neutrino transport. 
Since pairing quenches the rate by a large factor, all these 
nuggets freeze out with essentially the same baryon num- 
ber they had at formation. In general, the net result is 
that many nuggets survive with smaller masses, which 
could have not otherwise survived if paring had not op- 
erated. Therefore, any initial mass function of nuggets 
will be stretched towards the low-mass region after being 
partially evaporated. Note that this behavior is obtained 
for evaporating fluxes that may be many orders of mag- 
nitude larger than the very small values indicated by the 
chromoelectric flux tube models. 

We conclude that the survival of the nuggets (if 
formed) is much likely if they settle in a CFL state at a 
temperature T A = 57x (A/100 MeV) MeV, which may be 
true for the whole population. Thus, CFL prevents fur- 
ther evaporation/boiling and effectively freezes out the 
masses of the nuggets. A detailed numerical study of the 
whole evolution of the nuggets is desirable to address this 
issue. 



We have discussed in this work the effects of QCD pair- 
ing on the evaporation/boiling rates of quark nuggets 
assumed to be formed during the cosmological quark- 
hadron phase transition. These nuggets would be pro- 
duced at Tqcd ~ 150 MeV with maximum baryon num- 
bers A max ~ 10 49 (100MeV/T QC D) 2 corresponding to 
the horizon scale at that epoch. After formation, the 
nuggets are fragile because of the hot environment and 
may boil and/or evaporate into hadrons. The nuggets 
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